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Experimental data for central collisions of ^^®Xe on ^^®Sn at beam energies of (a) 32 MeV/nucleon, 
(b) 39 MeV/nucleon, (c) 45 MeV/nucleon and (d) 50 MeV/nucleon are compared with results 
calculated using a hybrid model. We use a transport model (BUU) to obtain the excitation energy 
per nucleon in the center of mass of the multifragmenting system. The canonical thermodynamic 
model is then used to determine the temperature which would lead to this excitation energy. With 
this temperature we use the canonical thermodynamic model to calculate various experimental 
data such as multiplicities of different composites, probability distribution of the largest cluster etc. 
Agreement with data establishes the validity of the model. 

PACS numbers: 25.70Mn, 25.70Pq 


I. INTRODUCTION 


Nuclear multifragmentation is an important phe¬ 
nomenon, the study of which can reveal reaction mech¬ 
anism in heavy ion collisions at intermediate and high 
energies [IHEl. Central collision fragmentation reactions 
around fermi energy domain are extensively used for pro¬ 
ducing neutron rich isotopes and for studying nuclear liq¬ 
uid gas phase transition. 

Different theoretical models have been already developed 
for throwing light on the nuclear multifragmentation re¬ 
action and for explaining relevant experimental data. 
The models are mainly classified into two categories: (i) 

S namical models p-Q and (ii) Statistical models [3- 
. The dynamical models are based on more micro¬ 
scopic calculations where the time evolution of projectile 
and target nucleons are studied. In statistical models 
the clusterization technique is nicely incorporated but 
the disadvantage of statistical model is that the calcu¬ 
lation are started by assuming some initial conditions 
(like temperature, excitation energy, freeze-out volume, 
fragmenting source size etc.). These condition are either 
parameterized or obtained from some experimental ob¬ 
servables. 

In this work we develop a hybrid model for explain¬ 
ing multifragmentation reaction around fermi energy do¬ 
main. We treat central collision only. Initially the exci¬ 
tation of the colliding system is calculated by using dy¬ 
namical Boltzmann-Uehling-Uhlenbeck (BUU) approach 
013 with proper consideration of pre-equilibrium emis¬ 
sion. Then the fragmentation of this excited sys¬ 
tem is calculated by Canonical Thermodynamical model 
(CTM) [^. The decay of excited fragments, which are 
produced in multifragmentation stage is calculated by an 
evaporation model [l^ based on Weisskopf theory. Dif¬ 
ferent observables like charge distribution, largest clus¬ 
ter distribution etc. are calculated by using this hy¬ 
brid model for ^^®Xe-|-^^®Sn reaction at different projec- 
tile energies and compared with experimental data (20| 


The idea of setting the initial conditions for a statistical 
model from a dynamical model is of course not new; see 
for example Barz et al [l^. In many Statistical Model 
of Multifragmentation (SMM), the initial conditions are 
fixed by some measured data. In our hybrid model the 
initial conditions for the thermodynamical model are set 
up almost entirely by the transport model calculation. 
The concept of temperature is quite familiar in heavy 
ion collision and it is a better observable (compared to 
energy) for studying liquid gas phase transition. One 
standard way of extracting temperatures is the Albergo 
formula M, where temperature is calculated from the 
measured isotopic yields (i.e. cold fragments). Another 
common technique for obtaining temperature is to mea¬ 
sure the kinetic energy spectra of emitted particles. But 
in this method, the effect of sequential decay from higher 
energy states, Fermi motion, pre-equilibrium emission 
etc complicate the scenario of temperature measurement. 
Our hybrid model calculation is an alternative method 
for deducing the freeze-out temperature and it bypasses 
all such problems. 


II. BASICS OF THE DYNAMICAL MODEL 


The hybrid model consist of three different stages: (i) 
Initial condition determination by dynamical BUU calcu¬ 
lation, (ii) fragmentation by canonical thermodynamical 
model and (iii) decay of excited fragments by evapora¬ 
tion model. 

We start our calculation when two nuclei in their respec¬ 
tive ground states approach each other with specified ve¬ 
locities. The mean field potential energy density is taken 
from Lenk-Pandharipande [2l[: 
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FIG. 1: (Color online) Evolution of test particles at (a) 0 fm/c, (b) 75 fm/c and (c) 200 fm/c in center of mass frame for 45 
MeV/nucleon ^^®Xe on ^^®Sn reaction. 


where first two term represents zero range Skyrme in¬ 
teraction and the derivative term does not affect nu¬ 
clear matter properties but in a finite system it pro¬ 
duces quite realistic diffuse surfaces and liquid drop bind¬ 
ing energies. This can be achived for A =-2230.0 MeV 
/m^, B=2577.85 MeV <7 =7/6, po = 0.16 and c=- 

6.5 MeV/m^/^ j2l|. We first construct Thomas-Fermi so¬ 
lutions for ground states [l^. The Thomas-Fermi phase 
space distribution will then be modeled by choosing test 
particles with appropriate positions and momenta using 
Monte Carlo. Each nucleon is represented by 100 test 
particles {Ntest = 100). We begin Boltzmann-Uehling- 
Uhlenbeck (BUU) model calculation to get the excita¬ 
tion of the fragmenting system. In the center of mass 
frame the test particles of the projectile and the target 
nuclei (in their Thomas-Fermi ground state) are boosted 
towards each other. The test particles move in a mean- 
field U{p{f)) (generated by the potential energy density 
of Eq.(l)) and will occasionally suffer two-body colli¬ 
sions when two of them pass close to each other and the 
collision is not blocked by Pauli principle. The mean- 
field propagation is done using the lattice Hamiltonian 
method which conserves energy and momentum very ac¬ 
curately [2l|. Two body collisions are calculated as in 
Appendix B of ref. except that pion channels are 
closed, as there will not be any pion production in this en¬ 
ergy region. Positions and momenta of the test particles 
are updated after each time steps (At) by the equations 


^ = -yrU{p{ri),t) 


dn 



* — 1 , 2 ,., {Ap + At)Ntest ( 2 ) 


III. EXCITATION ENERGY DETERMINATION 

We can calculate the excitation energy (E*) from pro¬ 
jectile beam energy {E{,eam) by direct kinematics by 
assuming that the projectile and the target fuse to¬ 
gether. In that case the excitation energy is E* = 
ApEbeam/iAp + At) where Ap and At are projectile and 
target masses respectively. This value is too high as a 
measure of the excitation energy of the system which 
multifragments. Pre-equilibrium particles which are not 
part of the multifragmenting system carry off a signifi¬ 
cant part of the energy. 



FIG. 2: Variation of energy of the central dense region (con¬ 
taining 80% of total test particles) with time obtained from 
dynamical BUU calculation for ^^®Xe on ^^®Sn reaction at 45 
MeV/nucleon. 
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FIG. 3: Left Panel indicates the variation of excitation energy 
per nucleon with projectile beam energy per nucleon obtained 
from dynamical BUU model. The Canonical Thermodynami¬ 
cal Model (CTM) can calculate average excitation energy per 
nucleon for a given freeze-out temperature, mass number and 
charge. Therefore to know the required freeze-out temper¬ 
ature corresponding to each excitation (obtained from BUU 
calculation) CTM is used. The variation of freeze-out tem¬ 
perature with projectile beam energy is shown in the right 
panel. 



FIC. 4: (Color online) Theoretical charge distribution 
(red dotted lines) for ^^®Xe on ^^®Sn reaction at (a) 32 
MeV/nucleon (b) 39 MeV/nucleon (c) 45 MeV/nucleon and 
(d) 50 MeV/nucleon. The experimental data are shown by 
black squares. 


To get a better measure of excitation of the fragment¬ 
ing system we need to do a BUU calculation where the 
pre-equilibrium particles can be identified and can be 
taken out to calculate excitation energy per nucleon. 
We exemplify our method with central collision reactions 
^^®Xe-|-^^®Sn at projectile beam energy 45 MeV/nucleon. 
Initially the center of ^^®Xe and ^^®Sn are kept at (lOOfm, 
lOOfm, 90fm) and (lOOfm, lOOfm, llOfm) respectively 
and and they are boosted towards each other along z 
direction. Fig. 1 shows the test particles at t=0 fm/c 
(when the nuclei are separate), 75 fm/c (the time when 
violent collisions occur) and 200 fm/c (almost all colli¬ 
sions are completed). From the figure it is clear that 
for t=200 fm/c some test particles are far distant from 
the central dense region. These fit the category of pre¬ 
equilibrium emission. In different multifragmentation ex¬ 
periments, it is observed that after pre-equilibrium emis¬ 
sion, 75% to 80% of the total mass creates the fragment¬ 
ing system [23l - l^ . Hence we choose the test particles 
which create 80% of the total mass (i.e. Aq = 198) from 
the most central dense region. Knowing the momenta of 
selected test particles the kinetic energy is calculated and 
from the positions of these selected test particles the po¬ 
tential energy is calculated by using Eq. 1. By adding ki¬ 
netic and potential energy the energy of the fragmenting 
system is obtained. Fig. 2 shows the variation of excited 
state energy of the central dense region (i.e. 80% of the 
total test particles) with time. Here total energy is always 
constant but as time progresses, pre-equilibrium particles 
having high kinetic energy, are escaping from the central 


dense region, therefore the energy of the central dense re¬ 
gion is decreasing. It is clear that after t=100 fm/c, the 
energy becomes independent of time. Hence, we can stop 
BUU calculation at any time after t=100 fm/c and con¬ 
sider the corresponding energy as excited state energy. 
To get the excitation we need to know the ground state 
of the fragmenting system. For this we use the Thomas 
Fermi method for a spherical nucleus of mass A = 198 
(80% of ^^®Xe-|-^^®Sn mass). Subtracting ground state 
energy from the calculated energy above the excitation 
energy is obtained. 


IV. COMPUTATIONS WITH THE 
STATISTICAL MODEL: EXTRACTION OF 
TEMPERATURE 

We have described above how from BUU we extract 
the mass, charge and the excitation energy of the frag¬ 
menting system. Our next task is to obtain the freeze- 
out temperature. The canonical thermodynamic model 
(CTM) 0 can be used to calculate average excitation 
per nucleon for a given temperature, charge and mass. 
Getting an excitation energy for a given temperature, 
mass and charge is described in detail in [^. We do the 
exploration for each beam energy. We will not repeat 
the formulae of CTM here but just mention that apart 
from neutrons and protons the following composites are 
included in CTM breakup. We include deuteron,triton, 
'^He and for heavier nuclei we include a ridge along 
the line of stability. The composites that follow from 
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CTM will further decay by evaporation. The details of 
how we do it can be found in ll^. 


V. RESULTS 

We have done calculations for the same ^^®Xe+^^®Sn 
pair for projectile beam energies 32, 39, 45 and 50 
MeV/nucleon. In each case, we have stopped the time 
evolution at t = 200 fm/c, and selected 80% of the total 
mass from central dense region for calculating the ex¬ 
cited state energy. Then subtracting the ground state 
energy the excitation is obtained. The variation of cal¬ 
culated excitation energy with projectile beam energy is 
shown in the left diagram of Fig. 3. From this excitation 
energy we find out the corresponding freeze-out temper¬ 
ature. Thus the freeze-out temperature for a given beam 
energy is obtained. This is plotted on the right side of 
Fig.3. 

To check the accuracy of our model, we have com- 
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FIG. 6: (Color online) Variation of average size of largest clus¬ 
ter with projectile beam energy obtained from hybrid model 
calculation (red dotted lines) for ^^®Xe on reaction. The 

experimental data are shown by black squares. 



FIG. 5: (Color online) Theoretical largest cluster probability 
distribution (red dotted lines) for ^^®Xe on ^'^^Sn reaction at 
(a) 32 MeV/nucleon (b) 39 MeV/nucleon (c) 45 MeV/nucleon 
and (d) 50 MeV/nucleon. The experimental data are shown 
by black squares. 

pared the theoretical results with experimental data. Fig. 
4 shows the comparison of charge distribution at pro¬ 
jectile beam energies 32, 39, 45 and 50 MeV/nucleon. 
With the increase of energy (i.e. increase of tempera¬ 
ture), fragmentation is more, therefore multiplicities of 
higher fragments gradually decrease. Fig. 5 represents 
the largest cluster probability distribution at different en¬ 
ergies. Since with the increase of energy breaking in¬ 
creases, the peak of the largest cluster probability distri¬ 
bution shifts towards the lower atomic number side and 
the width of the distribution gradually decreases. The 


variation of average charge of largest cluster {ZLargest) 
with projectile beam energy is shown in Fig. 6. In each 
case nice agreement between theoretical result and ex¬ 
perimental data is obtained. 


VI. DISCUSSION 

In this work we do a BUU calculation to get the exci¬ 
tation energy of the multifragmenting system produced 
in the central collision reactions around Fermi energy do¬ 
main, then do a CTM exploration to locate the temper¬ 
ature which will give this excitation. Once this temper¬ 
ature is fixed CTM is used to fit available experimental 
data. The agreement with data is pleasing. This work 
complements the work we did where we fitted the data 
obtained from the decay of projectile like fragments at en¬ 
ergies in the limiting fragmentation region. There [2^ we 
fitted the data using CTM with an assumed temperature 
profile first and later [l^ we showed that the tempera¬ 
ture profile is obtainable from BUU calculations. Our 
present calculations are not prohibitively computer in¬ 
tensive. One virtue of these calculations is that equation 
1 leads to reasonable values of binding energy of finite 
nuclei (even in Thomas-Fermi approximation) and real¬ 
istic diffuse surface without having to supplement the 
zero range Skyrme interaction with a finite range inter¬ 
action. Vlasov propagation for large nuclei when finite 
range interaction is present is very computer intensive. 
The other pleasing aspect is that the lattice Hamiltonian 
method |2l| gives remarkable accuracy in total energy 
and total momentum conservation in these calculations. 
For fragmenting system, we adopted the value of 80% 
of total mass from the experimental papers quoted in 
our paper. But our results (see Fig. 6) show that this 
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was a reasonable choice. We show a plot of {ZLargest) 
which agrees fairly well with data. Now {ZLargest) de¬ 
pends upon the size of the fragmenting system as well as 
the temperature of the fragmenting system. The larger 
the fragmenting system, the larger is the (ZLargest)- The 
higher the temperature, the smaller is {ZLargest)- Now 
the temperature also depends upon what percentage of 
nucleons are left out as pre-equilibrium particles. The 
value 80% we choose gives a combination of temperature 
and fragmenting mass that seems to be just about right. 
One could do a detailed best ’’fit” but this was not at¬ 
tempted. 

What we presented in this work did not involve any radial 
flow. One reason is that the collision energy being only 
about 50 MeV/nucleon, the initial compression is small 


so any radial flow must also be small. The best signa¬ 
ture for radial flow will be in the velocity distribution but 
we are only calculating multiplicity distribution. Neither 
CTM nor SMM can incorporate radial flow easily but in 
Lattice gas model, where flow is easily incorporated, it 
was found that even for significant radial flow, multiplic¬ 
ity distributions are hardly affected [^ . 
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